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Part I: 

Turbulence Modeling and Industrial 

Flows 


• Many industrial flows are turbulent; certainly in the markets 
that two of our codes, FLUENTand RAMPANT, are focused in. 

• Turbulence augments rates of mass, momentum and heat trans- 
fer, often by orders of magnitude. 

• Most combustion processes involve turbulence and often de- 
pend on it. 

• Choice of turbulence model dictates the accuracy of CFD pre- 
dictions. 

• There is still a large gap between the state-of-the-art and users' 
expectations and needs. 
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Challenges in Turbulence Modeling 


• Modeling the correlations: pu^u'j and pu-<£. 

- Closures based on the “eddy-viscosity” concept (indus- 
try’s most popular choice) 

- Closures based on transport equations (RSM) 

• Modeling an additional transport equation for a scalar quantity 
to fix the state of turbulence. 

- Most popular choice: the kinetic energy dissipation rate, 
e. 

- However, this equation is derived by continuum mechanics- 
based phenomenological considerations and intuition. 

• Modeling of the viscosity-affected, near- wall laminar sublayer. 

- Most popular choice: “Wall-functions” that bridge the 
turbulent field to the solid wall. 

- However, assumptions involved are not always right. 


Desirable Attributes of Turbulence 
Models in Commercial CFD Codes 


Accuracy and Universality 

— The range of applicability should be as broad as possible. 

— Applicable to complex geometries and unstructured meshes. 


• Economy 

— Mathematically simple. 

- Memory and CPU requirements should be moderate and 
affordable (model formulation and grid distribution re- 
quirements). 


• Robustness 

- Model should be able to solve a wide range of problems 
with little or no convergence problems. 

- Computationally efficient (fast execution speed and uses 
memory sparingly). 
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Turbulence Models in FLUENT 


• k-e model adequate for simple flows with no significant strain 
rates. 

• RNG k-e model for separated flows, flows with large streamline 
curvature, swirling flows, or flows with significant strain rates. 

• RSM recommended for swirling flows or highly anisotropic flows. 


k-e Model: 


Some Comments 


• Well-tested, used for over 20 years, limitations well understood. 

• It forms a good compromise between universality and economy 
of use for many engineering problems. 

• Subject to the inherent limitations of the Boussinesq’s hy- 
pothesis, i.e., isotropic eddy-viscosity and Newtonian closure 
(gradient-diffusion model). 

• Many assumptions are introduced in deriving the modeled equa- 
tions for the turbulent quantities, particularly the e-equation, 
making their fidelity limited. 

• The constants in the modeled equations are calibrated against 
simple benchmark experiments. 

• As a result, the k-e model performs poorly in flows with cur- 
vature, swirl, rotation, separated flows, low-Reynolds number 
flows, strongly anisotropic flows, etc. 
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Renormalization Group (RNG) 
Based k-e Model 


• Basic theory and derivation are described in Yakhot and Orszag 
(1986). Further details and applications are in Yakhot, Orszag, 
Thangam, Speziale, and Gatski (1992), Speziaie and Thangam 
(1992). 

• First introduced in a commercial code, FLUENT, in 1992. 

• The RNG method is essentially a scale-elimination technique 
that can be applicable to the Navier-Stokes and other scalar 
transport equations as well 

• Removal of successively large scales leads to differential trans- 
port equation models and associated formula for quantities such 
as the turbulent Prandtl/Schmidt number. 

• The basic form of the RNG-based k-e equations remains largely 
the same with the standard k-e model. But, the constants in 
the model equations are derived explicitly from theory. 

• The ^-equation ends up with an additional source term, a 
strain-dependent term. 

• The RNG model can be integrated directly to a solid wall with- 
out using ad hoc damping functions or damping terms used in 
many near- wall models. 


• High-Re form of the turbulence kinetic energy and dissipation 
rate equations derived by RNG procedure are: 
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P k = 2 vrSijSij is the kinetic energy production 
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RNG-Based k-e Model (Cont’d) 


• In the low- Re RNG model, a differential relationship exists 
between ^ and (Yakhot and Orszag, 1986). 

• The turbulent Prandtl/ Schmidt number is no longer a constant, 
and computed from relationships relating the local value of the 
number to the viscosity ratio (Yakhot and Orszag, 1986). 

• In these relations, as 0 — ► 1, a — ► a o (the low- Re limit) and 
as v — ► oo, a = a~ l — ► 0.7179 (the high- Re limit). Here: 

v — i / e ff/V 0 > where i/ eff - i/ 0 + i*r 
a = inverse turbulent Prandtl number (<r _1 ) 

Oq = inverse molecular Prandtl number (ctq 1 ) 


• In the low- Re regions, a* and <r e are obtained similarly from 
the Prandtl number relationships, with ao = 1.0. 

• The relationships ensure that in the high-Re number part of 
the flow where z> >> 1: 


i/ cff = i/j> = 0.085 — 

€ 

and the effective viscosity varies smoothly from the molecular 
viscosity to the turbulent viscosity. 

• The low-Re eddy-viscosity formula does not explicitly involve 
any geometric length scale, i.e., the distance from a solid wall 
used in the damping functions adopted by most low -Re near- 
wall models, which is a veiy convenient feature for calculations 
for complex three-dimensional geometries. 

• In collaboration with the originators of the RNG model, Drs. 
Yakhot and Orszag, the model has been extended to account 
for the effects of compressibility, swirl, rotation, and premixed 
combustion. 

• The RNG-based k-e model also works well with conventional 
and enhanced (non-equilibrium) wall functions available in Flu- 
ent Inc. codes. 
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The Reynolds- Stress Model in 
FLUENT 


• RSM solves transport equations for the Reynolds stresses: u-u'- 
(4 equations in 2D problems, 6 equations in 3D problems). 

• RSM is the level of modeling that has a well established track 
record of out-performing eddy- viscosity models in complex flows. 

• It is computationally more expensive and more inclined to di- 
vergence and stability problems. 

• The simple and most widely tested form of the Launder, Reece 
and Rodi (1975) form is used. 

• The interpolation technique for co-located grids of Rhie and 
Chow (1983) is used. 

• It offers the best choice for highly anisotropic flows. 


Example 1: 

Circle-to-Rectangle Transition Duct 


• Measured by Davis (1991). 

• Reo = 3.9 x 10 5 . 

• Solution Domain. 

- Upstream Inlet Boundary: x/D — —1.0 
— Downstream Exit Boundary: x/D — 8.0 
— A Quadrant of the duct modeled. 



STATION 1 
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Example 2: 
Cyclone Spar at or 


• Measured by Qing (1983). 

• RSM is used. 

• Cylindrical 55 x 23 x 41 grid. 
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Example 3: 

180° Bend of Square-Cross Section 

• Solution Domain 

- Upstream Boundary: 5.01 ? u from the start of the bend 

- Downstream Boundary: 5.0 D H from the end of the bend 

- A symmetric half of the duct modeled. 

• Mesh 

— Orthogonal 101 x 47 x 27 

- Distance from the wall « 0.01 Dm 
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Part II: 

Combustion Modeling 


• With environmental awareness. legislations on combustion- gen- 
erated pollutants such as NO x , SO x carbon monoxide, soot, 
unbumt hydrocarbons, etc. have become increasingly tougher. 

• Combustion simulation in industrial applications can help us to 
design combustors with higher efficiencies and lower pollutant 
emissions. 

• The combustion process involves some of the most complex 
phenomena such as chemistry, multiphase flow, turbulence, 
heat transfer and the interaction between these phenomena. 

• Here we focus on gaseous combustion in which the reactants 
may be mixed or non-mixed prior to flowing into the combustor. 
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Turbulence-Chemistry Interaction 


• Accurate simulation of turbulent combustion requires a thor- 
ough assessment of the way turbulence and chemistry interact. 
The reaction rate and flame structure primarily depend on this 
interaction. 

• In turbulent flames, chemical rates can be significantly differ- 
ent than those in laminar flames (sometimes several orders of 
magnitude), and the mean chemical rate is not the same as the 
rate calculated based on mean values of the various scalars: 

u02,...)*r(e u e 2 ,...) 

• Turbulent-chemistry interaction is best characterized by the 
Damkohler number which is the ratio of characteristic flow time 
to chemical reaction time: 


£>a = 


zz — a. 

Tr ~ if V 


• When Da << 1 chemical reactions are orders of magnitude 
slower than turbulent mixing and the influence of turbulence 
on reaction can be neglected. 

• When Da >> 1 chemical reactions are very fast and hence 
combustion is controlled by turbulent mixing. 


• At high Da we can exploit the laminar flame concept: turbu- 
lent flame is comprised of an array of laminar flames (flamelets) . 
Hence chemical rate expressions can be those obtained in lam- 
inar flames and the effect of turbulence can be characterized 
through the probability density function (pdf). 

• For turbulent diffusion flame, the pdf is usually expressed in 
terms of a scalar which can best characterize mixing, e.g., the 
mixture fraction. Since the rate of reaction is much higher 
than the mixing rate, we can assume that the reaction system 
is at equilibrium. The effect of turbulence is simply felt by the 
fluctuations in the mixture fraction. The mean value of any 
scalar in the flame is simply: 

o = Jo 

• For turbulent premixed flames the pdf is usually expressed 
in terms of a scalar which can best characterize the reaction 
progress, e.g., normalized temperature: 

r = /d r(c)p(c)dc 
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FLUENT Equilibrium Model 


• For turbulent diffusion flames we use a two-moment beta pdf 
and equilibrium data to calculate various thermo-chemical scalars 
in the flame. 

• To obtain equilibrium data we use the popular CHEMKIN li- 
brary of SANDIA, fully interfaced with our codes. CHEMKIN 
contains data on all important gaseous fuels, combustion inter- 
mediates and products as well as their properties. 

• We obtain the mean mixture fraction and its variance from 
their respective conservation equations: 



• To save computational time we calculate the integrals before 
the CFD calculations. 


Concluding Remarks 


• As of now, we provide our users with three turbulence models: 

— the “conventional” k-e model, 

— the ReNormalization Group model, 

— the Reynolds-Stress Model. 

• The Renormalization group k-e model has broadened the range 
of applicability of two-equation turbulence models. 

• The Reynolds-stress model has proved useful for strongly anisotropic 
flows such as those encountered in cyclones, swirlers and com- 
bustors. 

• Issues remain, such as near- wall closure, with all classes of mod- 
els. 

• Collaborative research with ICOMP will not only serve to fur- 
ther quantify applicability of turbulence models but may bring 
to market new ideas in the field of turbulence modeling for 
industrial flows. 
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